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Crossover from the chiral to the standard universahty classes in the conductance of a 
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The conductance of a quantum wire with off-diagonal disorder that preserves a sublattice sym- 
metry (the random hopping problem with chiral symmetry) is considered. Transport at the band 
center is anomalous relative to the standard problem of Anderson localization both in the diffu- 
sive and localized regimes. In the diffusive regime, there is no weak-localization correction to the 
conductance and universal conductance fluctuations are twice as large as in the standard cases. 
Exponential localization occurs only for an even number of transmission channels in which case 
the localization length does not depend on whether time-reversal and spin rotation symmetry are 
present or not. For an odd number of channels the conductance decays algebraically. Upon moving 
away from the band center transport characteristics undergo a crossover to those of the standard 
universality classes of Anderson localization. This crossover is calculated in the diffusive regime. 
Numerical simulations agree qualitatively with the theory. 
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I. INTRODUCTION 



Since the introduction of the scaling approach to the 
problem of Anderson localizationEa it is known that 
transport characteristics of a disordered metal are univer- 
sal, provided the disorder is sufficiently weak, the tem- 
perature sufficiently low so that quantum coherence is 
maintained over large distances, and the interaction be- 
tween electrons can be neglected. ppAn example is the 
phenomenon of weak localization,BQ a small deviation 
from Ohm's law for the conductance of a weakly disor- 
dered metal, which is suppressed by the application of a 
time-reversal symmetry breaking magnetic field. Though 
small, the weak-localization correction is universal in the 
sense that it does not depend on the shape of the sam- 
ple, nor on any other microscopic or macroscopic prop- 
erty other than its dimensionality and the presence or 
absence of time-reversal symmetry and spin-rotation in- 
variance. Another example is rtke phenomenon of uni- 
versal conductance fiuctuationsiffl The sample-to-sample 
fluctuations of the conductance of a disordered metal or 
semiconductor are of order e^jh with a prefactor that 
only depends on dimensionality and symmetry. Both 
the weak-localization correction and the universal con- 
ductance fluctuations are precursors of the true Ander- 
son localization, where as a result of destructive interfer- 
ence of multiple scattered quantum mechanical waves the 
dirty metal turns into an insulator for sufficiently strong 
disorder, or, in onje or two dimensions, for a sufficiently 
large sample size.Ll _ 

The original paper by Anderson)! and most of the 
effort devoted to the problem of Anderson localization 
since then, considers the case of a particle on a lattice 
with a random on-site potential (diagonal disorder) and 
non-random hopping amplitudes. In that case, one dis- 
tinguishes three universality classes, corresponding to the 



presence or absence of time-reversal and of spin-rotation 
symmetry. These three classes, are called orthogonal, 
unitary, and symplectic, respectively. Here, we will refer 
to these as the three "standard" universality classes. 

The electronic localization problem was soon general- 
ized to lattice models with randoruness in the hopping 
amplitudes (off-diagonal disorder). (This type of ran- 
domness was previously known from the description of 
phononsll^lO and narrow-gap semiconductors. t3) The lo- 
calization problem with off-diagonal disorder has received 
comparatively much less attentiau, although it has been 
known since the work of DysonE2l that random systems 
with off-diagonal disorder, but without diagonal disorder, 
can behave in a way dramatically different from that of 
systems wi1jli,|dia|gonal disorder only, or with both types 
of disorder .Il3ll3"ll3 For instance, the average density of 
states (DoS) for a one-dimensional chain with random 
nearest-neighbor hopping-sgas..faund to be singular at the 
center of thej-band, e — O.EjOlij According to the Thou- 
less formulajlj such a singular DoS implies that at ^r^A 
the conductance distribution be anomalous as wellJi3'EHl 
Wegner and Gade in Ref. (see also Refs. [6|j2|-|2|) 
found a two-dimensional counterpart to the singular be- 
havior of the average DoS within their analysis of a non- 
linear-tj model with a sublattice symmetry. Interest in 
the effect of off-diagonal disorder has revived in the 90's 
on two fronts. Motivated by quenched approximations 
to interacting theories such as the quantum Hall effect 
at half-filling or gauge approaches to high Tc supercon- 
ductivity, the random flux problem (a special case of 
off-diagonal disorder in which hopping amplitudes have 
a random phase only) has been extensively studied, al- 
though very 
has emerged 

vated by the close resemblance between the anomalies at 
zero energy induced by pure off-diagonal disorder in two 



.tie consensus on its localization properties 
A second thrust of activity has been moti- 
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dimensions and the nature of the plateau-transitions in 
the integer quantum Hall effect (IQHE):EII Both models 
might share the property that-all eigenstates are localized 
except at one special energy.cS 

The reason why the localization properties of the ran- 
dom hopping problem can depart from those of the stan- 
dard problem of Anderson localization is the existence of 
an additional sublattice symmetry in sjratcjms. with off- 
diagonal but without diagonal disorder OO^cSlj In that 
case, the lattice can be divided into two sublattices, such 
that the Hamiltonian changes sign under a transforma- 
tion where the wavefunction changes sign on one sublat- 
tice, but not on the other. As a result, the spectrum is 
symmetric with respect to a reflection about e = (i.e., 
eigenvalues appear in pairs ±e). The fact that the band 
center e = is a very special energy in the presence of 
the sublattice symmetry explains why anomalies in the 
DoS and the localization properties occur at precisely this 
value of the energy. When the energy moves away from 
zero, the effects of the sublattice symmetry on the spec- 
trum and the wavefunctions decreases and a crossover to 
the standard behavior takes place. The sublattice sym- 
metry is brokenJpy the presence of on-site disorder, long- 
range hoppingpJ or (in some cases) by periodic bound- 
ary conditions .123 Counterparts to this sublattice symme- 
try in other disordered systems or in quenched approxi- 
mation to interacting problems are numerous. They oc- 
cur in, the QCD Hamiltonian,Eilc3 random XY spin 
chains,E3 diffusion in randoro-environmentsj^ supersym- 
metric quaatum mechanics,0 non-Hermitean quantum 
mechanics,E3 and two-dimensional disordered models in 
the continuurruauch as Dirac fermions with random vec- 
tor potentials.EJ Following previous worlp..in this field, 
which adopted the nomenclature of QCD,e3 we will refer 
to the sublattice symmetry as chiral symmetry and will 
restrict our attention to random hopping problems with 
this symmetry. 

One-dimensional disordered systems with chiral sym- 
metry have been well-studied with all kinds of approaches 
and in various contexts (for references, see the previous 
paragraph) , and despite a continuing confusion about se- 
mantics, their localization properties can be considered 
well-understood. For two-dimensional systems the situa- 
tion is different (see Refs. and references therein). 
Reliable analytical and numerical results are notoriously 
hard to obtain, and no consensus has been reached to 
date, not even on some most elementary issues. In view 
of this controversy, it is particularly instructive to study 
the natural intermediate between one and two dimen- 
sions, the thick (or "quasi-one-dimensional") disordered 
wire. On the one hand, it shares the existence of both 
a localized and a diffusive regime of quantum transport 
with two dimensional disordered systems, while on the 
other hand, it allows for a controllable analytic treat- 
ment, just like the truly one-dimensional system. More- 
over, quasi-one-dimensional systems appear as a logical 
intermediate step in the finite-size scaling approa(±-|for 
numerical simulations in two and three dimensions.E3 



Localization properties at the band center of a quasi- 
one-dimensional quantum wire with off-diagonal disor- 
der were investigated in several previous publicatiai|is-.by. 
the authors, together with Simons and Altland.EjtH'EJ 
In those works we derived,a_chiral counterpart to the so- 
called DMPK equation,c3 l3 a Fokker-Planck equation 
that governs the distribution of the transmission eigen- 
values of a quantum wire without chiral symmetry. So- 
lution of the chiral DMPK equation for lengths beyond 
the localization length of the standard DMPK equation 
showed that there is no exponential localization if the 
number N of propagating channels is odd (including the 
one-dimensional case) , while the conductance decays ex- 
ponentially with length if N is even. This parity effect 
is strikingly similar to the sensitivity of the low-energy 
sector of a singjle. antiferromagnetic spin-A^/2 chain to 
the parity of iV,c3 on the one hand, or to the sensitivity 
of the low-energy sector of N coupiled antiferromagnetic 
spin- 1/2 chains to the parity of A^,c3 on the other hand. 
In the special case of the chiral Fokker-Planck equation 
without time reversal invariance (random phase quantum 
wire), it was possible to calculate exactly the crossover 
from the diffusive to the localized regime for all moments 
of the conductance and to verify the validity of the as- 
sumption of universality against a numerical simulation 
of the random flux problem.L3 The numerical simulations 
also confirmed that sufficiently far away from the center 
of the band, transport is governed by the standard uni- 
versality classes. 

A limitation of the approach relying on the Fokker- 
Planck equations for the transmission eigenvalues is that 
it cannot describe how the conductance distribution 
crosses over from the chiral to the standard universality 
class as e is tuned away from zero. In the renormalization 
group language, each Fokker-Planck equation describes a 
fixed point corresponding to a case of pure symmetry and 
the fixed points by themselves cannot be used to infer 
how the scaling flows take place between them. One pos- 
sibility to obtain information about the crossover energy 
and length scales below (above) which the physics is that 
of the chiral (standard) universality classes, is to study 
the DoS of a chiral quantum wireo However, unlike in 
the case of a one-dimensional wire, where the Thouless 
formula connects conductance and DoS, for a quasi-one- 
dimensional wire it is not possible to infer transport prop- 
erties from the DoS. In this paper, we use an alternative 
approach, developed by one of us for the study of trans=. 
mission through a random waveguide with absorption.cS 
Focusing on weak-localization corrections and universal 
conductance fluctuations, we compute how, in the diffu- 
sive regime, the conductance distribution of a quantum 
wire with random hopping crosses over from the chiral to 
the standard universality classes as the energy is tuned 
away from zero. We are not able to compute the crossover 
in the localized regime. Instead, for the localized regime, 
we consider the conductance distributions in the pure 
symmetry classes and compare to numerical simulations 
to establish the crossover scale and to verify the validity 
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of our predictions. 

The paper is organized as follows. In Sec. |l^, we define 
our microscopic model and derive the symmetries of the 
scattering matrix in the presence of the chiral s ymm etry. 
We then explain the scaling approach in Sec. III. The 
localized regime is studied in Sec. IV. Our main results 
are presented in Sec. 0, where we consider the crossover 
from the chiral to the standard universality classes in the 
diffusive regime. In Sec. we compare our theoretical 
predictions to a numerical simulation of a random hop- 



ping model on a square lattice. We conclude in Sec. VI] 



referred to as the orthogonal and unitary symmetry class 
and are labeled by the symmetry index f3 — 1 and 2, 
respectively. The case of broken spin-rotation symme- 
try with time-reversal symmetry is denoted /3 = 4 and is 
referred to as the symplectic class. When /3 = 4, the ele^. 
ments of the N x N matrix STm are real quaternions .l3 
The situation when both time-reversal symmetry and 
spin-rotation symmetry are broken reduces to the uni- 
tary class (/3 = 2) and will not be considered separately 
in this paper. We further assume that ST^ has a Gaus- 
sian distribution, with zero mean and with variance given 
by 



II. MICROSCOPIC MODEL AND SCATTERING 
MATRIX 

A convenient microscopic model that describes a single 
particle hopping randomly between two sublattices is de- 
fined. The symmetries obeyed by the scattering matrix 
associated to this Hamiltonian are derived. Assuming 
weak disorder, the microscopic model is approximated 
by a model defined in the continuum, for which the scat- 
tering matrix is explicitly constructed. 
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(2.3) 



A. Microscopic lattice model with chiral symmetry 



where 



In a general form, the Schrodinger equation for an iV- 
chain system with random hopping between two sublat- 
tices and without on-site randomness reads 

-e^{m)^T,^^{m + \)+Tl_^^{m-l). (2.1) 

For a spinless particle, ^'(m) is the A^-component wave- 
function where the index m labels the position along 
the chain. For a particle with spin- 1/2, "^{m) is the TV- 
component wavefunction made of spinors. In that case, 
the N X N hopping matrix Tm consists of quaternions .E3 
The system, and the allowed hopping matrix elements 
are depicted in Fig. 1(a). Note that the case of a square 
lattice with nearest- nei ghbo r random hopping is included 
in the general formula (2.1), see Fig. 1(b). 

In order to model transport, we consider a disordered 
region of finite length L — Ma, a being the lattice 
constant, and attach ideal leads with hopping matrix 
Tm = lAf on both ends, see Fig. 1(a). Following Refs. 

we draw the hopping matrices T,„ with m inside 
the disordered region from a distribution centered around 
the N X N unit matrix, 



Tm = exp(JT„ 



(2.2) 



We distinguish three symmetry classes depending on the 
presence or absence of time-reversal and spin-rotation 
symmetry. For a spinless particle (or for a spin-1/2 par- 
ticle in the presence of spin- rotation symmetry) , the hop- 
ping matrix 6Tm is real (complex) if time-reversal sym- 
metry is present (absent). These two cases are commonly 



j = PN + 2-(3- 



2(1-^) 
N '■ 



(2.4) 



and I is the mean free path. (Why t can be identified as 
the mean free path is explained below.) Here, the sym- 
bol f denotes the operation of complex conjugation for 
(3 = 1,2, whereas it denotes the opera±|ion of Hermitean 
conjugation for quaternions for /3 = 4.c3 We assume weak 
disorder, f ^ a. The parameter -q governs the relative 
randomness of the determinant of Tm- (See Ref. |4^ for 
the reason for its introduction.) 

We have chosen the statistical distribution ( p. 3D for 
technical convenience; it allows for an exact solution of 
the transport problem. As a justification for this choice, 
we recall that the transport properties do not depend on 
details of the microscopic model as long as disorder is 
weak, ^ 3> a, and the length L of the system is much 
larger than the mean free path. All properties of the mi- 
croscopic model are summarized in the two parameters £ 
and rj. [The proper value of the parameter rj depends on 
the details of the microscopic model under consideration. 
For instance, for the random flux modelEj (which is a spe- 
cial case of a random hopping modeU, rj = Q, while > 
in generic random hopping models .cd] To emphasize this 
universality, we compare our final results to numerical 
simulations for nearest-neighbor random hopping on a 
square lattice, cf. Fig. 1(b). 
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n = 3 



n = 2 



n = 1 



complex conjugate of any eigenfunction is itself an eigen- 
function with the same energy. (For P = 4,|-Gpmplex 
conjugation is meant in the quaternion sense.c2l) Since 
outgoing and incoming plane waves are interchanged un- 
der complex conjugation, we infer that time-reversal in- 
variance is represented by the additional constraint 



S*S=1 



2N- 



(2.8) 




n = 3 



n = 2 



n = 1 



m — 



m — 5 



FIG. 1. (a) Random hopping model as described by Eq. 
(2.1), for = 3. A disordered section of the wire (of length 
M = 4) is attached to ideal leads. Different chains are only 
coupled in the disordered region; there is no coupling between 
the chains in the perfect leads, (b) The nearest-neighbor ran- 
dom hopping model on a "square" lattice is a special case of 
the model considered under (a). 

In the leads on the left (L) and right (R), the 



Schrodinger equation (2.1) at energy e is solved by a 
sum of plane waves moving towards the disordered re- 
gion (denoted by a subscript i) and away from the sample 
(denoted by a subscript o) (see Fig. 2), 



m^ — ikma 



(2.5) 



Here < ^ < 7r/a, e — —2 cos ka, and and tpf^ (tp"^ 
and ■0°^) are A'^-components vectors containing the am- 
plitudes of the incoming (outgoing) plane waves in the 
left and right leads, respectively. The amplitudes of the 
ingoing and outgoing wav es are connected through the 
Schrodinger equation (2d) in the disordered region. This 
relation is formulated in terms of the 2N x 2N scattering 
matrix S^, 



Current conservation implies 

S^S = ^2N■ 



(2.6) 



(2.7) 



(Here and below we suppress the index e if only scattering 
matrices at the same energy are involved.) For the cases 
(3 — 1,4, i.e., if time-reversal symmetry is present, the 



L 



FIG. 2. Quantum wire with a disordered region of length 
L — Ma. Incoming plane waves are ^/j^'^i and ip^^i- Outgo- 
ing plane waves are ■4^°k and There are A'^ channels, 
i.e., n, n' = 1, ■ ■ ■ , A''. In a quasi-one-dimensional geometry, 
L = Ma > Na. 



The Schrodinger equation (2.1) has an additional sym- 
metry: the Hamiltonian changes sign under the transfor- 
mation 5'(m) (— )'"\1/(to). Correspondingly, for any 
realization of the disorder, the spectrum of energy eigen- 
values is symmetric about the band center e = 0. This 
symmetry, which originates from the fact that the dis- 
order preserves the bipartite structure of the lattice, is 
referred to as chiral symmetry. The chiral symmetry is 
a special attribute of random hopping between different 
sublattices; it is broken by e.g. on-site randomness or 
next-nearest-neighbor hopping. It is the chiral symme- 
try that is responsible for the anomalous transport prop- 
erties at the special e nersy f = of a quantum wire 
with random hopping.li3nJ'c2l To find the effect of the 
chiral symmetry on the scattering matrix, we note that 
the transformation '^(m) — s- (— )™^'(m) interchanges in- 
coming waves at energy s into outgoing waves at energy 
— e, and vice versa. Applied to Eq. (12. 61), this gives 



= s_ 



(The second equality follows from Eq. (|2.6|) at energy 
—£■) Taken together with flux conservation ( p.7| ), we thus 
find that the presence of the chiral symmetry results in 
the constraint 




(2.9) 



(2.10) 



for the scattering matrix S. Unlike the constrai nts o f 
flux conservation and time-reversal symmetry, Eq. ( 2.10| ) 
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involves scattering matrices at different energies. The ex- 
ception is the band center e = 0, where we find that S is 
Hermitian, 



(2.11) 



The scattering matrix is decomposed into four N x N 
subblocks r, r' and t, t', the reflection and transmission 
matrices, 



5 = 



t' 



t 



(2.12) 



The transmission and reflection matrices determine the 
transport properties of the quantum wire. They are re- 
lated to the conductance of the wire through the Lan- 
dauer formula, 



G 



2e' 



titH : 



2e2 
h 



and to the shot noise powerH 



P = 



■tr [th{i - th)] = 



■P, 



(2.13) 



(2.14) 



V being the applied voltage. (See Ref. ^ for more appli- 
cations to quantum transport.) A further decomposition 
of S follows from the polar decomposition of the matrices 
r, r' and t, t\ 

„_/V' 0\/ tanhX (coshX)-i\/V 
1^0 U J \{coshX)-^ -tanhX J \0 W 

(2.15) 

where W , V, and V arc N x N unitary matrices and 
X is an iV X A'^ diagonal matrix with real numbers Xj 
[j = 1, • • • , N) on the diagonal. In the presence of time- 
reversal symmetry, one has 



WW = V*V 



1 



N- 



(2.16) 



Chiral symmetry implies a relationship between the uni- 
tary matrices U, W , V, and V' at opposite energies. 



W 



't 



'T 

— £1 



(2.17) 



I n ter ms of the eigenvalues Xj, the equations (2.13) and 
( 2.14 ) for the conductance and the shot noise power read 



N 

E 



tanh^ : 



cosh 



E 



cosh 



(2.18) 



B. Continuum model with chiral symmetry 

For weak disorder (mean free path £ much larger than 
the lattice spacing a), we may replace the lattice model 



^JJ) by a continuum model. We linearize th e sp ectrum 
of the kinetic energy of Schrodinger equation (2J) in the 
close vicinity of the band center e = 0. Choosing a rep- 
resentation with left and right movers, we arrive at the 
continuum Schrodinger equation 

- eilj{y) = [^(73 (g) tNdy era (g) v{y) + a2® w{y)] ip{y). 

(2.19) 

Here ^ is a 2N component vector (elements of ip occur in 
pairs that correspond to left and right movers), v and w 
are N x N Hermitean matrices, and the cr^ (/Lt = 1,2,3) 
are the Pauli matrices. In the presence of time-reversal 
symmetry w (v) is (anti)symmetric. The continuum limit 
has been taken along the chains only; discreteness is 
maintained in the transverse direction through the N 
components of ip. The Fermi velocity has been set to 
one. The randomness in the hopping amplitudes has 
been translated to the matrices v and by means of 
the identiflcations 



ST,„ - ST 



t 

m+l 



6T - ST 



m+l 



h.C. 
h.C. 



w{y). 



(2.20) 



With the choice (2_^), the disorder in v is statistically 
independent from the disorder in w. Both v and w are 
Gaussian distributed with zero mean and with variance 



{wi.j{y)wki{y')'') 



_ P5{y-y') 

2(/3-l)(l-r?) 

I3N 
P6{y-y') 



(2.21a) 



2(1-77) 



^^k^Jl + —Jf-^a^jk 



(2.21b) 



The symmetries (flux conservation, time-reversal, and 
chiral symmetry) of the scattering matrix in the contin- 
uum model are the same as for the lattice model. (Note 
that in the continuum model, the chiral transformation is 
represented hy ip aitp. The chiral symmetry then fol- 
lows from the fact that ai anticommutes with the Hamil- 
tonian.) 



III. SCALING APPROACH 

The idea0 behind the scaling approach to the the- 
ory of localization in a quantum wire is to calculate how 
the scattering matrix S of the quantum wire changes if 
a thin slice is added to the disordered region [see Fig. 
3(a)]. Here we are mostly interested in the eigenvalues of 
the matrix product tH — 1 — r'^r, i.e., in the parameters 



Xj of the decomposition ( p. 15 ). Hence, it is sufflcient to 
consider the reflection matrix r, and calculate how it is 
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changed upon the addition of a thin sHce. This change 
follows from the composition law 



t[{l 



(3.1) 



that gives the reflection matrix of two scatterers 1 and 
2 in series, in terms of the reflection matrix of the right 
scatterer (2) and all reflection and transmission matrices 
of the left scatterer (1), see Fig. 3(b). 



(a) 



SL 




L 



(b) 



























^ '"2 



FIG. 3. (a) A thin slice of length SL with a < (5L < ^ < L 
is added to the left of the disordered region of length L. (b) 
Two disordered regions 1 and 2 with scattering matrices Si 
and S2, respectively, in a quantum wire. 

If applied to a quantum wire, the only input in this ap- 
proach is the statistical distribution of the transmission 
and reflection matrices ti, t'l, ri, and r{ of the thin slice. 
The width SL of the slice is taken much smaller than the 
mean free path £, so that the change of r is small as well, 
although SL must remain large compared to the lattice 
spacing a for the continuum limit to be a good approxi- 
mation. Then, the scattering matrix Si of the thin slice 
can be calculated in the second-order Born approxima- 
tion from the Schrodinger Equation (2.19). The result is 



-W + ^[V,W], 
ti= 1 4 
r[= W 



iV- 



ieSL, 



where 



V 



iV ^ 



SL 



2 ^ 



^W^ + ieSL, 



(3.2a) 
(3.2b) 
(3.2c) 
(3. 2d) 



dyv{y), 



W 



SL 



dyw{y). 



Here we neglected terms that are of order (SL)'^. 
also ignored the y-ordering of the integrals in Eq. 



as it does not affect the statistical distribution of 5*1 in 
view of the delta- function co rrela tion of the random po- 
tentials V and w.] Using Eq. (2.21) for the distribution of 
the random potentials v and w, we find that the matrices 
V and W are Gaussian distributed with zero average and 
with variance proportional to the width SL of the thin 
shce, 



PSL 



2(/3-l)(l-,7) 



PSL 

2(1-^7) 



(3N 



(3.3a) 



(3.3b) 



Equations ( |3.lD -( |3.3| ) define the sc aling approach. They 
are exact for the c ontin uum model (2.19) with the statis- 
tical distribution ( 2.21 ) of the random potentials, which 
in turn was d erived from the random hopping lattice 
model (2T,2^) in the limit of weak disorder. A differ- 
ent c hoic e for the distribution of the hopping matrices in 
Eq. (2.3) would have led to different statistical proper- 
ties of the scattering matrix for a thin slice. However, as 
we will verify in Sec. ^ by numerical simulations, such 
differences are irrelevant in the sense of the renormaliza- 
tion group, i.e., they disappear for sufficiently long wires 
(longer than the mean free path £). 

Note that the reflection probability N~^tr r\r^ of a 
thin slice has average 



N-\trrlr^) = SL/£, 



(3.4) 



which justifies our choice that i is the mean free path. 

In terms of the matrices V and W, upon addition of 
the thin slice, the reflection matrix r changes according 
to 



r r + Sr, 



(3.5a) 



with 



5r = 2ieSLr — W + rWr — i{Vr — rV) + rWrWr 

- \ {W^r + rW^ + V^r + rV^) + VrV. (3.5b) 

We have not included terms of order VW as their con- 
tributions vanish upon disorder averaging. 

Several observations can be made already on the level 
of the evolution equation ( 3^ ) , in combination with the 
Gaussian distribution (3.3) of the matrices V and W. 
First, the distribution of r is symmetric under a change 
of sign, r ^ —r. This implies that the average of any odd 
function of r must be zero, for all values of the energy e. 

Second, at the band center e = 0, th e chir al symmetry 
implies that r is Hermitian, cf. Eq. (2.11). The Hcr- 



miticity is broken by the first term in Eq. (3.5b), which 
is proportional to the energy. 
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Third, the distribution of r is invariant under transfor- 
mations r ^ U rll'^ , where U is an orthogonal (unitary) 
N X N matrix for (3 = 1 (2). For zero energy, where r 
is Hermitian, this implies that the distribution of r de- 



pends on its eigenvalues tanhxj only, of. Eq. (2.15). As 



was shown in Refs. in this case, the scaling flow 

can be represented in terms of a Fokker-Planck equation 
for the distribution P{xi, • • • , xn] L), 



N 



dL -fi ^ dx, V'^ N I dx, 

l,j=l ^ / J 



J = I sinh(a;; - Xk)f. 



(3.6) 



fe<; 



( |3.6[ ) and (^j) to compare the transport properties for 
£ = and e far away from 0. (A comparison for the case 
of broken time-reversal symmetry only has already been 
given in Ref. ^.) In t he d iffusive regime we start from 
the evolution equation ( |3.5| ) directly, in order to include 
the £-dependence of the transport properties. Knowledge 
of the crossover as a function of e will allow us to specify 
what is meant by "e sufficiently far away from 0" , and 
hence when the standard DMPK equation (3.7) replaces 
the special Fokker-Planck equation (3.(:) in the random 
hopping problem. 



IV. LOCALIZED REGIME 



Away from the center of the band, r is no longer Her- 
mitian, and its distribution depends on both eigenvalues 
and eigenvectors. However, for s sufficiently far away 
from (this notion will be made precise below), the chi- 
ral symmetry has no effect on the scattering matrix, and 
P{xi, ■ ■ ■ ,xn; L) obeys the Fokker-Planck equation for 
the standard orthogonal, symplectic, or unitary symme- 
try classes, the so-caJJpi Dorokhov-Mello-Pereyra-Kumar 
(DMPK) equation,y'E 



dP 
'dL 



1 



N 



2{f3N + 2 - /?)£ ^ dxj 



dxi 



J = JJ^ I sinh 2xk \ JJ^ | sinh^ xi — sinh^ Xk \ 



(3.7) 



k<l 



There is no parameter i] in the DMPK equation; the 
presence of the parameter rj is special for the case of 
chiral symmetry at the band center s — 0. In the lan- 



guage of the Fokker-Planck equation (3.6), rj controls the 
relative strength of the diffusion of the center of mass 
X = {xi -I- • • • -I- xn)/N compared to that of the relative 
coordinates Xj — x. 

The most impo rtan t difference between the Fokker- 
Planck equations (3.6) and (|3.7[ ) are the symmetries of 
the Jacobians J. In Eq. ( |3.6| ), i.e., at the band cen- 
ter £ = 0, J is invariant under a simultaneous transla- 
tion Xj —^ Xj + Sx and under a simultaneous reflection 
for all j. [The translation invariance decouples 

{xi + --- + xn)/N 
and hence calls for 



the motion of the "center of mass" x 

from the relative coordinates Xj — x, 

the presence of the para met er r/ in Eq. ( |3.6| ) .] In the stan- 
dard DMPK equation (3/7), i.e., for energies e far away 
from the band center, J is invariant under a reflection 
Xj ~* —Xj for each j separately; there is no longer trans- 
lation invariance. It is the absence of this "local" reflec- 
tion symmetry at £ = that is responsible for anomalies 
in transport properties at £ = 0. In the remainder of this 
paper, we describe these in more detail, focusing on the 
distribution of the conductance in the localized regime 
L ^ N£ and on the quantum interference corrections to 
the conductance in the diffusive regime £ <^ L <^ N£. For 
the localized regime, we use the Fokker-Planck equations 



Differences between the conductance distribution at 
the band center £ = and away from £ = are most 
pronounced in the localized regime L 3> N£. Away from 
the band center, the conductance decreases exponentially 
with length, as is the case in the standard orthogonal, 
symplectic, and unitary classes. At the band center, how- 
ever, the exponential decrease of the conductance is only 
observed if the number of channels is even, while for an 
odd number p£ channels the conductance decreases only 
algebraically.EJ 

Exact calculations for the moments of the conductance 
in thc-statidard symmetry classes have been obtained for 
all /3,E3~E3 while for the chiral s ymm etry classes governed 
by the Fokker-Planck equatian. (3_^) only exact results for 
f3 = 2 and 77 = 1 are known E3 While we do not know of 
a way to extend our exact analysis of Ref. |3^ to the cases 
of orthogonal and symplectic symmetries, it is still pos- 
sible to extract the conductance distribution deep inside 
the localized re gim e L ^ N£ using the approximation 
scheme of Refs. p3|j57| , |5^ . This is done here. We are thus 
able to compare the average and variance of the conduc- 
tance, and the average and variance of its logarithm at 
and away from the band center £ = for the orthogonal, 
symplectic, and unitary symmetry classes for all values 
of rj. 

Our starting point is the Fokker-Planck equation (p.(^), 
which we rewrite in the form 

9P 1 l-rj\ d fdP ^^dfl 



N ) dxi \dxj dxj 



with the "potential" 

N N 

= — In I sinh(a;fe — x/)]. 



(4.1) 



(4.2) 



fc=l Z=fc+1 



Equation (4.1) has the interpretation that as L increases. 



fictitious particles with the coordinates Xj perform a 
Brownian motion subject to the repulsive two-body po- 
tential 17. Since f2 has a hard core, we may assume that 
xi < X2 < ■ ■ ■ < xn iov all L. In fact, as a result of their 
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repulsive interaction, the distances between the Xj^s will 
grow with increasing length, until eventually for suffi- 
ciently large L 



and 



a;i ^ X2 <C • • • <C xn- 
Then we may approximate 

ax-i 



(4.3) 



(4.4) 



and find that Eq. ( [4.l[) is solved by a Gaussian distribu- 
tion for the Xj, 



N 



P{xi, ■ ■ ■,xn;L)(x exp <! ^ (x^ - J 




(4.5) 



Here, the N x N matrix has the entries {EM)ij = 
1, and the channel-dependent "localization length" 
reads 



6- 



7£ 



PiN + 1^2j)- 



(4.6) 



For comparison, in the standard orthogonal and uni- 
tary symmetry classes, the probability distribution 
P{xi, ■ ■ ■ , Xn', L) in the localized regime is also given 
by a Gaussian of the type (4.5), but with 77 = 1, 7 = 
2{PN + 2- (3), and = {(3N + 2~ + /3j - p)H 

In the localized regime L ^ Ni only the Xj th at ar e 
closest to contribute to the conductance, cf. Eq. ( ^.18 ). 
For even N, they are xn/2 a-nd X(isi/2)+ii both of which 
are an average distance 



{xN/2) = -(•T(JV/2) + l) 



L 
1' 



J' 



(4.7) 



away from zero. The length scale f serves as the localiza- 
tion length for even N. For odd N, the conductance is 
determined by only one eigenvalue, a:(7v+i)/2, which has 
zero average. 



(a;(Ar+l)/2) = 0. 



(4.8) 



The presence of the eigenvalue X( jv+i) /2 with zero average 
is responsible for the absence of exponential localization 
in this case. 

The average and variance of the conductance and the 
average and variance of its logarithm follow from the 
probability distribution (O) . For even TV the results are, 
with an accuracy ©(L^/^^for the logarithms displayed, 




(In ,9) 



var In g 



The latter result shows that, in the localized regime, the 
conductance distribution is well approximated by a log- 
normal distribution; unlike the average conductance g 
itself, which has fluctuations that are much bigger than 
the average, its logarithm In g provides a good character- 
istic of the ensemble. 

For odd N, there is no exponential localization. The 
conductance has a broad distribution, which is neither 
characterized by the (average of the) conductance nor its 
logarithm. 



Pig) 



exp 



^(l_i_!i) ^arccoshV 



1-') 



(4.10) 



With this distribution and up to corrections of order 
L°/^°, the average conductance decays algebraically, 



1 — 1] 

N 



{g')=o(g)' 



(4.11a) 
(4.11b) 



while the average of its logarithm grows proportional to 
L^/^ rather than L, 



var liig= 



1 - 



1 - 



l-r]\ L 



N 



(4.11c) 
(4.11d) 



Away from the band center e = 0, the conductance 
distribution follows from the standacd DMRK equation 
It is close to log- normal, witho^tSal^ 



L 

In var g = ln{g), 
2L 

4L 



ln(.g> 

1 var g 
(ln.g) = 



In 



var \ng 



6- 



(4.12a) 
(4.12b) 
(4.12c) 

(4.12d) 



Hg) 



f3 



N 



L 1 (L 



In var g = ln(5r). 



up to an accuracy of 0{Lr /S^). Here the localization 
(4.9a) length for the standard symmetry classes is given by 

(4.9b) U^iPN + 2- 13)1. (4.13) 
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The most striking difference in the conductance dis- 
tribution appears for odd N, where the absence of ex- 
ponential locahzation at e = is contrasted with the 
exponential decay of the conductance for e ^ 0. How- 
ever, also for an even number of channels, there is an 
important difference. At e = 0, the localizatio n le ngth 
£ K N£ is /3-independent for large N, cf. Eqs. (2.4) and 
(4.7), while the localization length ^st ~ PN£ away from 



t he b and center is proportional to /3 for e 7^ 0, cf. Eq. 
(4.13). Hence, upon moving away from the band center, 
the localization length increases by a fact or (3 . [The mean 
free path does not depend on e, see Eq. (3.4).] 

The absence of a /3-dependence for the localization 
length at the band center may be related to the anomaly 
in the DoS for random hopping models at that energy. In 
Ref. |4^, it was shown that in the absence of time-reversal 
symmetry the DoS p{e) near zero energy has a pseudo- 
gap, p(e) cx e| lne|, while in the presence of time-reversal 
symmetry p has a logarithmic divergence, p(e) cx jlnej. 
We conclude that, upon breaking time-reversal symme- 
try, the decrease in the DoS available for transport, can- 
cels the suppression of destructive interference respon- 
sible for the increase of the localization length in the 
standard case. 

The average and variance of the conductance in the 
localized regime are dominated by rare events, where the 
smallest Xj is close to zero (corresponding to a transmis- 
sion coefficient close to unity). For wires without chiral 
symmetry, approxi mati on of P{xi, . . . , xn] L) by a Gaus- 
sian similar to Eq. ( [1.5| ) fails for Xj close to zero because 
it does not account for the repulsion between Xj and its 
mirror image —Xj, cf. Eq. ( |3.7[ ). While it does not affect 
the leading 0{L) behavior of h\{g) and Invar 5, this fail- 
ure shows up in the subleading logarithmic terms in Eqs. 
( [4.124[l.l2b| ) which are different from what one would 
have obtained from a Gaussian distribution for the x, . 



[The results quoted in Eqs. (4.12a 4.12b) above. 
from an exact solution of the DMPK equation! 
In the presence of the chiral symmetry, there is no re- 
pulsion between Xj and — x, , so that the approximation 



(4.5) remains valid for Xj close to zero. In t his respect, 
we remark that the logarithmic terms in Eq. (^^), which 
were obtained with the help of Eq. (4.5) indeed agree 
with the exact solution of Ref. ^ for the case (3 = 2. 



V. DIFFUSIVE REGIME 



In the diffusive regime £ ^ L ^ N£, the effects of 
quantum interference do not take such a dramatic form as 
in the localized regime. The typical conductance of any 
sample is given by the classical Ohm's law, g = N£/L, 
and does not know of quantum mechanical phase coher- 
ence, the presence or absence of time-reversal symmetry, 
or, as we shall see below, the presence or absence of chiral 
symmetry. The role of quantum mechanics, and hence 
the role of the symmetries of the microscopic Hamilto- 



nian in this regime is confined to small corrections to the 
average conductance and to its sample-to-sample fluctu- 
ations. In spite of their smallness, these corrections are 
of prime importance, as they are a universal signature 
of quantum phase coherence, their size being determined 
by the fundamental symmetries of the system only. They 
do not depend on microscopic properties of the quantum 
wire, nor on its macroscopic characteristics, such as mean 
free path, width, or length. 

The two corrections are referred to as "weak- 
localization" and "universal conductance fluctuations" . 
The former is a small correction Sg (Sp) to the ensem- 
ble averaged (dimensionless) conductance (5) (shot-noise 
(p)) that is suppressed if time-reversal symmetry is bro- 
ken by a magnetic field. For a standard quantum wire, 
it reads 



6g = 



(3-2 
3/3 



Sp = 



(3- 



45/3 



(5.1) 



Since it signals the first departure from Ohm's law, the 
weak localization correction to the conductance is pre- 
cursor to the exponential suppression of the conductance 
in the localized regime. The universal conductance fluc- 
tuations refer to the sample-to-sample fluctuations of the 
conductance, which have variance, 



vai g 



2 

15^' 



(5.2) 



The breaking of time-reversal (spin-rotation) symmetry 
reduces the conductance fluctuations by a universal fac- 
tor of V2 (2). 

In this section we calculate those quantum corrections 
for the case of a quantum wire with random hopping. 
Our calculations are inspired by the approach of Mello 
and Stone ,E3 who have derived and solved scaling equa- 
tions for the moments of the conductance in the stan- 
dard universality classes from the DMPK equation in the 
limit of large N. We consider both the quantum correc- 
tions for the pure symmetry classes, corresponding to the 
Fokker-Planck equations ( |3.6D and (3.7) at e = and e 
far away from 0, respectively, and for the intermediate 
regime, where the crossover between the two symmetry 
classes takes place. Since in the latter case no Fokker- 
Planck equation for the transmission eigenvalues Xj is 
available, a modification of the approach of Ref. ^ is 
needed, which is based on the more funda men tal scal- 
ing equation for the reflection matrix r, Eq. (3.5), rather 
than on a Fokker-Planck equation for the transmission 
eigenvalues pXp . Such a method was proposed by one of 
the authorsLj in the context of the transmission through 
a random waveguide with absorption. Bel ow, we adapt 
this method to the present case (Sec. V A ), and present 
solutions for th e ch iral symmetry classes at the band cen- 
ter £ = (Sec. VB) and for the crossover from the chiral 
symmetry classes to the standard universality classe s as 
e moves away from the band center e = (Sec. VC). 
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A. Scaling equations 

Although we are primarily interested in the statistics 
of the transmission matrix t, and in particular in the (di- 
mensionlcss) conductance g — trtH and shot noise power 
p = tr — tH)\ , we find it more convenient to formu- 

late our scaling equations in terms of the reflection matrix 
r. Once we know r, unitary of the scattering matrix al- 
lows us to find the transmission properties without much 
effort. 

Before we write down the most general scaling equation 



J 



for a trace of an arbitrary product of the reflection ma- 
trix r and its Hermitian conjugate, we would like to focus 
on the scaling equation for tr r^r in order to demonstrate 
the method and the approximations involved. Addition 
of a thin disordered slice to a disordered wire causes a 
smal l cha nge r ^ r + Sr to the reflection matrix r, see 
Eq. (3.5b). Hence, upon addition of this slice, the t race 
trrV changes to tr {r + 6ry {r + Sr) . Using Eq. ( |3.5lj ) for 
6r, we thus find, up to 0{5L) [Recall that the variance 
of W is of order SL, so keeping terms up to 0{SL) means 
up to 0{W^)], 



6 tr rV = — tr [r(l — r^r)W + W{1 — r''r)r^ 

-(1 - rrt)VF(l - r^r)W + r(l - r''r)WrW + Wr^W^(l - rV)rtl 



(5.3) 



All terms that involve the disorder potential V in Eq. ( 3.5b ) canceled due to the cyclicity of the trace. Next we 
perform a disorder average over W and over the reflection matrix r of the wire of length L. We thus find 

= ([tr (1 - rV)]^) - (trrtrr(l - rV)) - (trr^ trr^(l - rr^)) 

p oL 

• ^ ~ '^(tr (1 - rV)(l - rV - - r^^)) - '^^^—^{tr (1 - rr^)(l - rV - - r^^)). (5.4) 



/3 



below Eq. (3.5). The initial condition at L = corre- 
sponds to perfect transmission, i.e., (trrV) = 0. The 
solution is easily found. 



(tr rV) 



Ns 



-0(7V"), 



(5.6) 



where s = L/£. This solution corresponds to Ohm's law 
for the conductance g ^ N — tr r^r, 



(.9) 



N 



I 

Finally, we take the limit 5L < C ^, and replace the finite Here we have used the fact that the average of the trace of 
differences on the l.h.s. of Eq. (5.4) by differentials. an odd product of r's and r^'s is zero, see our discussion 

It is apparent that the scaling equation obeyed by 
(trrV) is not closed: On the r.h.s. traces and products of 
traces of up to four reflection matrices appear. Closure 
requires an infinite family of scaling equations, and can- 
not be achieved on the level of scaling equations for the 
moments, but only with the help of the Fokker-Planck 
equation for the transmission eigenvalues Xj in the cases 
of pure symmetry. However, for lengths L ^ N£ it is 
possible to decouple this infinite set, and to find a solu- 
tion order by order in L/ {N£). Formally, this decoupling 
scheme proceeds along the lines of a large- expansion: 
In addition to the explicit factors A^ in Eq. (5.4), each 
trace contributes a factor A^. Further, we assume that, 
to leading order in A^, the average of a product of traces 
equals the product of the averages. As we will see be- 
low, corrections correspond to a (co) variance of traces, 
and are of order A^*^. Similarly, if we have a product of 
n traces, we can expand in cumulants, where an nth cu- 
mulant will turn out to be of relative size A^^"". Such a 
decoupling scheme is kno"wpi,to work for the case of the 
standard DMPK equation)^ and its consistency can be 
verified from the scaling equations for traces and prod- 
ucts of traces that we derive in this section. 

Let us now see how the scaling equation for (tr r^r) 
decouples in this large-A^ decoupling scheme. Recalling 
that 7 is of order A^, cf. Eq. (2.4), we thus find that the 
r.h.s. of Eq. (5.4) is of order A^i.e., 

idL (trrV) ^ N - 2(trrV) + ^ (tr rV)% 0(A^") 



0{N°) 



(5.7) 



To this order in A^, the result is entirely classical. The 
average (trrV) (and hence {g)) does not depend on the 
energy e nor on the presence or absence of time-reversal 
symmetry. The dependence on time-reversal symme- 
try shows up through the term proportional to on 



the r.h.s. of Eq. ( |5.4| ), which is of order A^. It is this 
term in the scaling equation that gives rise to the weak- 
localization correction to the conductance. The scaling 
equation for (trrV) does not contain an explicit energy 
dependence. Instead, the energy-dependence shows up 
through the appearance of the traces like trr^ or trrV'^ 



in Eq. (5.4) in the weak-localization correction. Such 



(5.5) 



traces that contain different numbers of r's and r^'s 
strongly depend on energy, as can be seen from the scal- 
ing equation of, e.g., (trr^). 
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^dLitrr^) - ([tr (1 - r^)]^) - 2(trrtrr(l - r^)) + l!£I^(trr2 



J 

2-P 2(1-77) 



(tr(l-r2)(l-3r2)). (5.8) 



/3 (3N 

With the same decoupHng scheme as before, we find a closed scahng equation for (trr^) up to 0{N), 

edLitir^) = TV - 2(1 - 2ie£){tYr^} + ^(trr^)^ + 0{N°), (5.9) 

which has the solution 

(trr^) = /V |l - 2iee + 2^e£{i + e£) cot 2v/rf(i + e£)s | \o{N°). (5.10) 

One verifies that for e — > 0, the average (trr^) equals the average (tr r^r) that we computed above, since for zero 
energy one has r — r'' . One also verifies that for e£ ^ 1 the average (trr^) approaches zero, as is the case in the 
standard symmetry classes. 

We are now ready to discuss the scaling equations for the trace of the product of an arbitrary number of reflection 
matrices and for the product of such traces. Hereto we write ro = r and ri = , and define 

=trrj, ...rj„, (5.11) 

where the indices jk can take the values or 1. We define the symbol R without indices as R = N. We also define 
products of traces through the symbols 



(i) (i) 

where rij denotes the n-tuple , . . . , z„y . 

Proceeding along the same lines as above, we then find that the scaling equation for a single trace is given by (see 
Fig. I) 



jdL{Rn-u) ^ 



V //? R 2(1-7?) 

/ J \-'^jk---jl-'^jr--jnjl---jk ' p ^jk---jdk---jljn---jl ^3k---jljl---jnjl---jk 



l<k<l<n 



V //? R ^^R 2(1-7?) 

/ , \ ^Jfc + l---ii-l^ji + l---jnil---jfc-l "T ^jk + l---jl-ljk-l---jljn---jl + l O TVf ^jk+l---jl-ljl + l---jnjl---jk-l 

2-/3o 2(1-7?) 



- \^ / R R I _f}li_lliR 

/ , \^''jk---jl-l-^^jl + l---jnjl---jk "T p ^jk---jl-ljk---jljn---jl + l ^3k---jl-ljl + l---jnjl---jk 

- V //? R +^/? 2(1-7?) \ 

/ ^ \-'^jk + l---jl-'^jvjnjl---jk-l^ n ^jk + l---jljk-l---jljn---jl ajvr ^jk + l---jljl---jnjl---jk-l / ■ {O.-LO) 

l<k<l<n' ^ PI 

Here, it is understood that jn+i = Ji, Jo = 3n- Moreover for 7i > Z = fc + 1 > 1, Rj^^^...j^_^ = trlAr = N ^ 
Rjk+i---ji-ijk-i---jijn---ji+i = Rjk-i---jijn---jk+2^ ^^'^ ^jfc+ 1 ■ ■ ■ ji - 1 ji+ 1 ■ ■ ■ ji ■ ■ ■ - 1 = "^jfc +2 ' ' ' ji ' ' ' jfc - 1 ' rcspectivcly, whcrcas 
when k — 1 and Z = tt, Rji^i---j„ji---jk-i = trl^v = N, Rjk+i---ji-ijk-i---jijn---ji+i = Rjk+i---ji-iji+i---jnji---jk-i = 
Rj2---jn-n respectively. Note that there is a one-to-one correspondence between contributions involving a product 
of two traces, say, 

Rjk---ji-Rjvjnji---jk = ti' {"Tj^ ■ • ■ rji) tr (rj, • • • Tj^rj^ ■ ■ - rj^.) , 
and contributions arising in the presence of time reversal symmetry 

'p ^ik---3i3k---iii^i---3i — 'p [{"""jk ■ ■ '^iJ (^jfc ■ ' ■'^ji'^jn ■ ' ■'^ii)] = tr (rj^. • • • J (rj, • • -Tj^Tj^ ■ ■ -Vj^) 

or due to the randomness in the determinant of the hopping matrices 

2(1-7?) 



13N 



'-^3k---jljl---jnjl---jk — [('"jfc ■ ■ ■ fjl ) {fjl ■ • • 1^jn''^jl ' ' ' ''jfc) 
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For products of traces, we find 



7^ 
J 



Ql {Qm- 



E 



fni---n3_inj + i---n, 



E \'3ni---nfc_infc+i---n,„ini+i---n„-f'nfcni / 1 (5.14a) 



l<k<l<m 



where ( ji/ l3)dLRnj is given by the r.h.s. of Eq. (|5.13| ) with omission of the angular brackets for the disorder averaging, 
and where F^^ ^ J2T=iYll=i fk,i with 



2(1 -r;) 



(3N 



+ 7? 2(1 -ry) 

"'"-'^Jfc+l---lmIl---ifc_ljl + l---j„Jl---Jl_l ^ p ^lk+l---lm.ll---lk-ljl-l---jlin---jl + l ^ «fc + 1 ' ' ' im »1 ' ' ' Jfc - 1 ^Jl + 1 • ■ ■ J 1 ' ' ' Jl - 1 

-7? -^R I 2(1-^7) 



2-/? 



2(1-77) 



^ik+l---imil---ik-ljl---jnjl---jl p ^ik+l---im.il---ik-ljr--jljn---jl f^]\j ^ik + l---imil---ik-l-^jl---jn.jl---jr 



(iN 



(5.14b) 



r 



Here we denoted m. — ii, . . . ,im and n — ji, . . . , jn- 

Below, we are interested in averages and (co)variances 
of traces of an even number of reflection matrices up to 
order N'^. In both cases, the terms proportional to (1 — 77) 
do not play a role. For the ave rage of a single trace, this 
is immediately clear from Eq. (5.13). To see this for the 
( co) va riance of two traces, some further inspection of Eq. 
( ^.14 ) is needed. First, 77 appears explicitly in the quan- 



t ity fm m multiplying a product of two traces, see Eq. 
( ^.14b| ). A priori, the leading contribution, which is ob- 
tained by replacement of those traces by their averages, 
i s of t he same order [0{N)] as the other terms in Eq. 
( 5.14b| ). However, as ttt, and n are even, each of the two 
traces multiplying (1 — 77) contains an odd number of re- 
flection matrices, so that their averages vanish. Hence, to 
leading order in N, the contribution from the term pro- 
portional to (1 — 77) vanishes. Second, 77 appears implic - 
itly through the derivative {'ji/ /3)dLRn- in Eq. ( 5.14a ). 
Again, to leading order in N, its contribution vanishes, 
and one is left with a term of relative size 7V~^. 



It should be mentioned that Eqs. ( |5.13D and ( |5.14| ) 
can be extended to the case in which a weak stagger- 
ing of the hopping am plitude is present in the micro- 



scopic model, cf. Eqs. (2.1 



2.4). (How to generalize the 



Fokker- Planck equation (3.6) to include dimerization was 
shown in Ref. see also Ref. ^.) Weak staggering of 
the hopping amplitude is implemented by requiring that 
the disorder potential W has the Gaussian distribution 
with variance (3.3b) and mean (Wjk) = ^^A6jk- Here 
A measures the strength of the dimerization along the 
chain direction. With weak dimerization, Eq. (5.13), say, 
is modified by the addition on the r.h.s. of the contribu- 
tion AX;Li(^j»= -j. -jfc - Rjk+i-j^-]k^i)- We see that 
the scaling equations now couple traces over an even and 
odd number of reflection matrices as is expected since the 
probability distribution of W is not anymore symmetric 



about W ^0,ci. Eq. (|3.5bD . 



Equations (5.13) and (5.14) are the central results of 
this section. These equat ions are more general than the 
Fokker-Planck equations (3.6) and (3.7) in the sense that 
they are valid both at the center of the band e = and 
in its proximity. Their limitation is that they can only 
be solved in the diffusive regime L <C Ni. In particular, 
they cannot be used to probe the localized regime (in 
contrast to their counterparts in the problem of a wave 
guide with absorption, see Ref. ^). The next two sub- 
sections are devoted to a solution in the diffusive regime. 
The ca se of pure chiral symmetry (e = 0) is considered 
in Sec. VB; t he en ergy dependence of the solution is dis- 



cussed in Sec. VC 
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FIG. 4. Diagrammatic representation of Eq. (5.13). Each 
circle corresponds to a reflection matrix r or rT To cal- 
culate the increment of {Rji---j„) (single box containing n 
open circles) one chooses a pair of (filled) circles. As in- 
dicated in (b) the same circle can be chosen twice. Over- 
lapping or nested boxes represent multiplication of traces. 



a2J3J4-"'J4J5JlJ2 / 1 



Thus, {Riij203iiK) ^ {RiiRjzjij&jiJ^K 
{RjiR^Ji)^ (-Rj2J3-Rj5JiJ2>> and (-Rj3j4-Rj4J5Ji>. ^rc repre- 
sented by (a), (b), (b'), (c), (d), and (e), respectively. 



B. Diffusive regime in the chiral limit 



The general scaling equations ( 5.13 ) and ( ^.14 ) sim- 
plify considerably at the band center £ = 0. At the band 
center, the scattering matrix is Hermitian, and hence 
r = r^. Restricting our attention to single traces and 
products of two traces, we find the scaling equations 



7^ 
1 



dL{trr" 



n7 

T 

n 

^2 



trr" 



ji-i 



k=0 

(3 2(1-,,) 



^(trr"-'=+i trr'^+i) + ^Y1 {^{"^^ r'''-^^^ ii r^-^) - 2 (tr r""''' tr r*^) 



fc=i 



7^ 



^^(trr^trr") = 



7^ 
1 



(3 (3N 
(tr (9Ltr r'")r") + (trr™(9Ltrr")) 
2mn{l — if) 



(n+ 1) (trr"+2) + (n - 1) (trr^^^^ - 2(n - 1) (trr") 



217171 



(3N 



(trr™-i(l - r^) trr"-i(l - r^)) . 



+ ^(trr"+"-2(i_^2)2^ 



(5.15) 



(5.16) 



Here {j£/(3)dLtT:r^ is the r.h.s. of Eq. (|5.15D with the 
omission of the diso rder averaging brackets. (If n — I, 
the last term in Eq. (5.15) should be omitted.) [Alterna- 
tively, one could have used the Fokker-Planck equation 
( |3.6[ ) to derive these scaling equations. Both methods 
agree as we have verified explicitly.] 

The average and variance of the conductance g = 
N — tr can be c ompu ted by straightforward solution 



of E qs. KWj ) and ( |5T6| ) using the decoupling scheme of 
Sec. VA. The result is, up to corrections of order N^^, 



(g) = 

var g = 



N 
s + 1 

4 
15^ 



2-/3 



13 (l + s)3' 
6s + 1 



(S + 1)6 



(5.17a) 
(5.17b) 



where, as before s — L/£. For the derivation of these 
results we needed the following intermediate results. 



(trr-*) = 
(trr^) = 



Ns^{3s^ + 8s + 6) 
3(l + s)4 ' 



177s2 + 180s + 75) 



15(1 



up to corrections of order N'^ and 

1 



(tr r tr r) 
(trr tr r^) 



2 

3^ 
2 

15/3 



1 



5s-' 



15s2 24s + 4 



(s + 1)6 



up to corrections of order N^^. 

In the diffusive regime £ L N£ we observe that 
the variance of the conductance at e = is twice the 
value taken in the standard case, for s far away from 



r 



0, cf. Eq. (5.2). The result that the presence of the ex- 
tra chiral symmetry leads to a doubling of conductance 
fluctuations was found previously for the random flux 
model (corceaponding to our case f3 = 2) from numerical 
simulationg23 and fro rtuan exact solution of the Fokker- 
Planck equation (3.6).EJ The factor two decrease in the 
fluctuations as the chiral symmetry is broken, is rem- 
iniscent of the factor two decrease of the conductance 
fluctuations upon breaking time-rejiepjal symmetryca or 
upon breaking a spatial symmetry.Ello 



According to Eq. ( 5.17a ), application of a magnetic 
field has an effect on the average conductance, but this 
effect vanishes in the diffusive limit £ <^ L <^ N£, i.e., 
s 3> 1. In other words, there is no weak- localization 
correction to the conductance in the diffusive regime. 
It is instructive to note a coincidence between the /3- 
dependence of the average conductance (g) and the /3- 
dependence of the localization length ^ which was con- 
sidered in the previous section. In the case of the ran- 
dom hopping model at zero energy, there is no weak- 
localization correction, and the localization length ^ does 
not depend on the presence or absence of time-reversal 
symmetry. On the other hand, without chiral symme- 
try (for large energies), the negative correction to the 
average conductance ior (3 — 1 foreshadows the localiza- 
tion transition, which occurs on a length scale ^st that is 
proportional to /3, i.e., localization takes place twice as 
fast without than with a time-reversal symmetry break- 
ing magnetic field. Absence of weak-localization correc- 
tion to the conductance had been pointed out by Gade 
and Wegner in their study of a (two-dimensional)-, non- 
linear-CT model implementing the chiral symmetry.Ell (See 
also Rcf. 



Finally, notice that there is no precursor in Eqs. (5.17) 
of the even-odd effect seen in the localized regime. This 
agrees with the exact solution for = 2, where it was 
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found that the even-odd effect is non-perturbative in the 
expansion parameter L/7V££3 

To order TV, the average conductance is the same as in 
the case of a wire without chiral symmetry. Differences 
show up only to order N'^, where we find that there are 
no weak-localization corrections for the chiral case. This 
is not a coincidence that is limited to the average of the 



J 



2-/5 



n 



conductance g = trtH. It extends to the averages of 
traces of arbitrary powers of r or t. To see this and in 
order to allow for a more detailed comparison to the case 
where chi ral sy mmetry is absent, we rephrase the scaling 
equation (5.15) in terms of the transmission matrix t. In 
the limit of large N, one thus obtains 



Wl (tr (th)") =-nY^ (ti- (i^i)"^^"") (tr (i^^)"') + " ^ (tr (t^i)"-™) (tr (t^t)™) 



rn— 1 

n 



m — l 



i2n + 1)^ (tr (ttt)"+i) + 2n2i-i^ (tr [t^r) + ©(TV"). 



(5.18) 



r 



For a quantum wire without the chiral symmetry, the 
leading 0{N) contribution to (t r (tH)") is precisely the 
same as the first line of Eq. (5.181). The weak- localization 
correction proportional to — differs in the standard 
case from the second line of Eq. ( |5.18| ) as_it reads 
_n2 2^(tr(tti)») + n{n - l)2^(tr (ttt)»+i).i| Hence, 



whereas the solution of Eq. (5.18) is the same as for an 



ordinary quantum wire to leading order in N, 

N 



{tritHr)^-B{n,l/2) 



(5.19) 



where 
function)^ 
2n(tr (th) 



] = T{x)r{y)/T{x + y) is the beta 
the combination (2n -I- l )(tr { thT+'^) - 
on the second line of Eq. ( ^.18| ) conspires 
with the coefficient Bin, 1/2) to insure the disappear- 
ance of a weak-localization correction for all averages 
(tr {tHY) in the presence of the chiral symmetry. As 
a corollary, we find that the average density of the trans- 
mission eigenvalues Xj 



p{x) 




Xj) 



(5.20) 



has no weak-localization correction as well. 

With these results, it is little work to compute the 
average shot noise power {p) = {ivtH{l - th)) and its 
weak-localization correction, 



«4 



1 



1 



+ 1 (s+l)^ 



+ 



2-/3 

/3 



Is' 



(5.21) 



_3(l + s)3 3(l-hs)6_ 

Just like in the case of the conductance there is no weak- 
localization correction in the diffusive regime £ ^ L <ti 
Nl. 



long length L of the quantum wire, its transmission prop- 
erties will flow to those of the standard symmetry class. 
This flow is governed by a crossover length scale so 
that for L <^ £g, the transmission properties are still alike 
those in the chiral symmetry class, while for L ^ they 
resemble those of the standard symmetry class. We dis- 
tinguish three possible regimes where this crossover can 
take place: 

• The crossover takes place in the ballistic regime, 

4 

• The crossover takes place in the diffusive regime, 

^ < 4 < NC, or 

• The crossover takes place in the localized regime, 
4 ^ Ni. This regime cannot be treated with the 
methods used in the paper. For the case iV = 1 
of a single-channel quantum wire, this regime has 
been studied in Refs. |T3|,f5llT8|j65|j20|. 



The full set of scaling equations ( ^.13| ) and ( 5.14 ) can be 
used to describe the flrst two regimes (and the intermedi- 
ate region between them). Although the solution of the 
scaling equations is straightforward — within the large- 
N decoupling scheme, the scaling equations are linear 
ordinary differential equations that can be solved one-by- 
one (see appendix ^) — it is a quite cumbersome task, 
and m any expressions get quite lengthy. [The expression 
( |3.10| ) for (trr^) is the only example whose solution can 
be represented by a one-line equation.] To simplify our 
presentation and to save the reader from those lengthy 
expressions, we focus on the regime <C -^e <C N£, where 
the crossover takes place inside the regime of diffusive 
dynamics. 

The len gth scale for the crossover can be identified 
from Eq. ( ^.IC ) as 



Crossover between the chiral and standard 
universaHty classes 



4 = 



(5.22) 



For any nonzero energy e, the chiral symmetry of Eq. 
(2.11) is broken. Hence one expects that for a sufficiently 



(Here we have neglected e£ with respect to i in y/i + e£. 
This is consistent with our focus on the regime ^ £.) 
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We note that e is none but the Thouless energy for a £ <^ we then find that the solution of the scahng equa- 

diffusive process with diffusion constant Vf£ (having mo- tions takes a relatively simple form. For the average and 

mentarily reinstated the Fermi velocity Vf) in a system variance of the conductance g we find up to 0{N^) 
of linear size fg. Using the hierarchy of length scales 



(.9) 
var g 



N£ 
2 



2-/3 



1 z coth(z*(T) + z* coth(z(7) 
3 



2 

15/3 ^ ]3 



+ c.c. 



Here we defined z = 



1 + i and a = L/ 
Nl 2-/3 r 1 



3cr4 



/3 



3zcrcoth(z(T) — 2 i 

16^^ + 8r72 sinh2(z*(7) 

For the average of the shot-noise power we find up to 0{N^) 
(3z - 2z*ct2) coth(z(T) 



45 



24cr3 



4a-2 sinh^(z*(T) 



(5.23) 
(5.24) 

(5.25) 



For the derivation of these results, we needed the following intermediate results, all up to corrections of order A^", 



(tr (rtr)3) - 
(trr^) = iV- 
(trrV^) = N - 
{trr'^^r^) = N - 
(tr(rV2)2) = N - 

(tr [r^rfr"^) ^ N- 



Ne 23 
~ 15ct' 

N£ 



4 
N£ 

Ne 



-z* coth(z*(T), 



(4(t2 + z)z* coth(z*cr) 



4a^ 



2crsinh^(z*cr) 



4(7 — z* coth(zcr) — zcoth(z*o') 



2(7^ 

(16cr^z* -|-8cr^z-f z*)coth(z*cr) zcoth(z(T) Scr^ - 2z(T coth(z*cr) + 5i 



lecr'i 4a4 8cr3 sinh^ (z*a) 

(gOcr^z* 40zct2 -f 3z*) coth(z*cr) 40(7^ - 6zcr coth(z*cr) + 3i' 



48cr3 sinh^(z*cr) 



(5.26) 



In the limit a ^ 0, corresponding to L <C (e, the 
weak-localization corrections to the conductance and the 
shot noise power and the conductance fluctuations ap- 
proach their values for the chiral symmetry class cf. Eqs. 
( ^.17 ) and (5.21). For cr » 1, corresponding to L :$> 
one verifies that the values corresponding to the standard 
symmetry class are recovered. 

In this subsection, we have described the effect of a 
finite energy by a crossover length scale fg. For the con- 
ductance and the shot noise power the limits of large 
and small e correspond to the limits of L large or small 
cor npare d to £/,. However, upon inspection of Eq. ( 5.10| ) 
or (5.26|) one observes that traces like trr^ that contain 



of L. 

We are now ready to define what is meant by "e suffi- 
ciently large" in the crossover from the chiral symmetry 
class to the standard symmetry class. As far as quantum 
interference corrections to the transmission properties are 
concerned, the results of this subsection show that "e 
sufficiently large" corresponds to the inequality of length 
scales L ^ or equivalently, e ^ £/L'^. However for 
reflection traces like tr , a much more strict criteria is 
needed, 



< ^, or £ < 



In the next section, these 



different numbers of r's and r^'s, do notr-approach their 
large-energy limits (trr^) = as L — > oo.E3 The origin of 
this difference is that the reflection matrix is dominated 
by (interference of) paths that only enter a distance of the 
order of a mean free path i into the quantum wire, while 
the conductance and the shot noise power depend on 
quantum interference throughout the entire wire. Hence, 
as long as £e ^ £, the finite energy cannot alter the in- 
terference of most paths that contribute to r. Hence, to 
judge whether the finite energy is relevant for the traces 
of refiection matrices, one has to compare £s to £ instead 



criteria, as well as the functional forms ( 5.23 ), ( 5.24 ) for 
the crossover will be compared to numerical simulations. 



VI. NUMERICAL SIMULATIONS 

In this section we report on numerical simulations of 
the conductance of a quantum wire with random hor 
only, and compare them with the theory of sections [ 

The simulations are for the random hopping model on 
a square lattice, described by the Schrodinger equation. 



-tm-l,]-\\^rn-l,j " *m,i; || V'm+l.j , 



(6.1) 
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where 'ipmj is the wave function at the lattice site {m,j). 
A site is labeled by the chain index j = 1, . . . ,N and by 
the column index m. We impose open boundary condi- 
tions in the transverse direction, im,o;± = tm,N:± ~ 0. 
The system consists of a disordered region (0 < m < L), 
coupled to the left and right to perfect leads (m < 1 and 
m > L). In the leads, the longitudinal and transverse 
hopping amplitudes are t^j-w = 1 and tm,j-i_ — t, where 
< t < 1. With this choice, there is a window of ener- 
gies — l + t<e<l — t around the band center, where 
the number of transmission channels does not depend on 
energy (and equals the number of chains N). In the dis- 
ordered region, the hopping amplitudes are taken from 



a distribution centered around the values t„ 



,r,\\ 



= 1 and 



tm,j;± = t for the leads. We consider two types of ran- 
domness, that we refer to as the real random hopping 
and random flux models. 

• In the real random hopping (RRH) model, the hop- 
ping amplitudes t,n,j;± and tm,j-\\ are chosen uni- 



formly and independently in the intervals —t{l 



S) < t„ij;± < t{l + S) and 1 — S < < 



1 + 6, respectively, where S measures the disorder 
strength. A uniform magnetic field with a flux 0pi 
through each plaquette is modeled by multiplica- 
tion of im,j;|| with a Peierls phase e^^^'^p^^^~^\ 

In the random flux (RF) model, the longitudinal 



hopping amplitudes t^, 



,r,\\ 



— 1, while the trans- 



verse hopping amplitudes tm,j;± are complex num- 
bers tmj;± = ie**"' 
such that the fluxes 



Here the 9,^ j are chosen 



are inde- 



pendently and uniformly distributed in the interval 
—TTp < 4>m,j < Trp, where p is a measure for the 
strength of the disorder. 

In the random flux model, the parameter 77 = 0, see Ref. 
PI ; in the real random hopping model the precise value 
of 77 is not known. However, nonzero rj (of order by 
assumption) will only give rise to corrections of relative 
order 1/iV, for the average conductance), which 

can be neglected for large N. Since the statistics of the 
conductance in the RF model in a quasi-one-dimensional 
geometry has been studied extensively in Ref. ^ at and 
away from the band center e = 0, we restrict our atten- 
tion here to the crossover as a function of energy. 

Th e wavefunctions that solve the Schrodinger equation 
(3.1) at energy e can be written as 



^ 1 

= E [e''=""sin(g.j>^L(^) 



e-'^^-^sm 



('ZA'+i-.J^f (^)] (6.2) 



in the left lead and as 



sm 



+ e"=""sin(g,j)'Af (^)] 



in the right lead, where the wave number fc^ > is 
determined from e = — 2cosfc^ — 2tcosgi/ with = 
'iTv/{N + 1). With this parameterization, the definition 
of the scattering matrix and its symmetries are the 
same as in Sec. O. 

For each realization of the disorder, the dimcnsionless 
conductance g is computed from the Land pjcr formula 
( pl3| ). The recursive Green's functionEjSE3 method is 
used to calculate S^- (Application of the method to the 
random hopping or random flux models is discussed in 
Ref. |3^.) Our numerical simulations use the parameters 
t = 0.6, 5 = 0.2, and p = 0.3, for the RRH and RF 
models. 



A. Localized regime in the RRH model 



In the localized regime, the even-odd effect manifests 
itself most dramatically. Taking an average over 2 x 10''^ 
realizations of the disorder, we have computed the mean 
and variance of g at the band center e = for the RRH 
model with = 20 and 21, and with and without a 
time-reversal breaking magnetic field, see Fig. |. The 
magnetic field corresponds to a flux 0pi = 8 x 10^'' per 
plaquette, or 1 flux quantum per 50 lattice spacings 
along the chain, so that time-reversal symmetry is broken 
for all but the shortest wire lengths shown in Fig. ||. For 
odd N{= 21) both {g) and varg decrease algebraically 
whereas they decay exponentially for even A^(= 20). We 
observe that, for odd (even) N and flxed L, {g) and var g 
are larger (smaller) in the presence of a magnetic field, 
13 = 2, tha n without, 13 — 1, in agreement with Eqs. (4.9) 
and (4.11). Note that for small L, varg is L-independent 
for the chiral unitary class, while y&v g decreases linearly 
with L for small L in the chiral orthogonal class. Similar 
L-dependencies for small L have been obtained for the 
standard symmetry classes, see Ref. pil 
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FIG. 5. Mean (a) and variance (b) of the conductance for 
N — 20 and 21 at the band center £ = with and without 
magnetic field in the RRH model. Averaging over 2 x 10* 
realizations of disorder is performed. 



FIG. 6. Localization length ^ as functions of energy e 
for (^pi = 0, 2 X 10"*, 4 X 10"*, and 6 x 10"* in the RRH 
model with A'^ = 20 chains. Averaging over 500 realizations 
of disorder is performed. 



Results for the crossover from the chiral universality 
classes to the standard ones as a function of energy are 
shown in Figs. ^ and|^. Figure ^ shows the energy depen- 
dence of the lo calization length ^ = — 21inii_>oo L/{\ng) 
[cf. Eq. ( [4.9c| )] for N = 20; Fig. ^ shows numerical 
data for the ratio C = — limi_>oo(lng)/var Ing. Here, 
the averages were taken over 500-1000 realizations of 
the disorder and magnetic fields corresponding to fluxes 
(j>pi = 2, 4, 6 X 10"* per plaquette, respectively, have been 
used. 



In the absence of a magnetic field, ^(e) shows non- 
monotonic behavior with a maximum around s ~ b x 
10^^, while, within 10%, the localization length ^ is 
the same in the chiral orthogonal class (e = 0) and in 
the standard orthogonal class (e > 10"* for the choice 
of parameters in the simulations), in agreement with 
As we discussed in Sec. the fact that 
0) — £,{e ^ 0) in the absence of a magnetic 



Sec. |IV. 

as 

field, could be interpreted as the result of a cancella- 
tion of two effects: The presence of an extra symmetry 
at the band center (the chiral symmetry), which tends 
to make ^ shorter than away from the band center, and 
the enhancement of the_DoS at the band center, which 
tends to make ^ larger.c3 Apparently, these two com- 
peting effects are not balanced in the crossover region, 
thereby giving rise to the non-monotonic energy depen- 
dence to ^ displayed in Fig. ^. Such a non-monotonicity 
of i^(e) is reminiscent of the non-monotonic voltage de- 
pendence of differential conductanqe-|found in a normal- 



metal / superconductor microbridge .1 



In the presence of a magnetic field, ^ increases by a fac- 
tor ~ 1.8 in the crossover from the chiral- unitary to the 
standard unitary symmetry class, which is slightly less 
than the factor 2 predicted in Eq. (4.13). Note that the 
increase in localization length is most rapid around the 
same energy scale Sc ~ 5 x 10^^ for which reaches its 
maximum, ^c, in the absence of a magnetic field. More- 
over, Ec is related to by Thouless relation ^ \J ^l^c 
where the mean free path i is obtained by dividing the 
localization length ^(e = 10"*°) by TV = 20 in Fig. |. 

A plot of the ratio C — — limi^oo(ln5)/var In g is 
shown in Fig. ^. In the standard symmetry classes, C 
takes the univjeusal value C — while in the chiral 
classes one haald 



C = 



/3N/2 



N+{1- 



0{1/N). 



2r?) 



(6.4) 



The data shown in Fig. ^ confirm that C — 0.5 in the 
standard symmetry classes (corresponding to e > 10~* 
for the parameters of our simulation). However, for the 
chiral symmetry classes, a 20% discrepancy with Eq. 
(U) is found. 
While the numerical simulations for the localized 



regime qualitatively confirm the theory of Sec. tV, quan- 
titative agreement is only up to ^ 20%. As a possible 
source of this discrepancy, we point to the fact that the 
simulations are done for an appreciable disorder strength 
d = 0.2, while the theory is derived for weak disorder, 
corresponding to (5 ^ 0. Hence, the system cannot be 
considered truly (quasi) one-dimensional, and corrections 
from two-dimensional dynamics on shorter length scales 
need to be taken into account. Another cause of the 
observed discrepancies could be the uncertainty of the 
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precise value of rj for the RRH model. While we believe 
that rj should not affect the conductance distribution sig- 
nificantly for large N, it remains difficult to make a quan- 
titative assessment of finite- iV corrections as long as 77 is 
unknown. At this moment, we are not aware of a direct 
way to obtain ry from the numerical simulations. 
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FIG. 7. Ratio C — — limL_oo(ln(?)/varlngi versus e for 
(f)pi = and 4 x 10""* in the RRH model with A*' = 20 chains 
and L = 4 X 10^. The disorder average is taken over 500 
realizations. 
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B. Diffusive regime in the RRH and RF models 



FIG. 8. Variance of g for the RRH model versus the total 
flux (j) through the disordered wire (a) and versus energy e 
(b). In (a), varg is shown for two values of the energy, £ = 
and e = 0.005, corresponding to the presence and absence of 
chiral symmetry; in (b), varg is shown for two values of the 
flux (j) = Q and </> = 8, corresponding to the presence and 
absence of time-rever sal sy mmetry. The dashed line in (b) is 
the theoretical curve ( |5.24| ), with cr^ = 1.5 x 10*e. The simu- 
lations are performed for A'^ = 45, L — 800, and an averaging 
over 10* reahzations of the disorder is performed. 



We next consider the crossover from the chiral sym- 
metry classes to the standard symmetry classes in the 
diffusive regime I <^ L <^ In Figs. |[ and || we show 
numerical simulations of the average and variance of the 
conductance as a function of the energy e and the mag- 
netic flux (j) = L{N — l)(/>pi through the disordered part 
of the wire. The simulations are performed with A'^ = 45 
and L = 800, in order to ensure that the conditions 
£ <C i <C ^ for diffusive transport and N <^ L ior quasi- 
one-dimensionality are both met. The ensemble average 
is taken over lO'' samples. 



Numerical results for the variance of the conductance 
in the RRH model versus (j) and e are shown in Fig. ^. The 
numerical data of vax g versus (j) (Fig. |^) agree within 
10% with the theoretical predictions va.r g = 4/15 (2/15) 
for (/) = and var^ = 2/15 (1/15) for </> > 1 for e = 
(e ^ 0). The crossover between the orthogonal and uni- 
tary classes happens for 0^1, both with chiral symme- 
try (e = 0), and without (e = 0.005). The e-dependence 
of var g is sh own in Fig. ||b, together with the theoretical 
result ( 5.24 ), where we fitted the crossover energy scal e 
that enters into the definition of cr = L/£^, cf. Eq. (5.22). 
Again we find quantitative agreement well within 10%. 
(The fact that the numerical data for (p = are below the 
theoretical curve can probably be attributed to the sup- 
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pression of var 5 as L approaches the locaHzation length 
^, see the remark in the discussion of Fig. ||.) 
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FIG. 9. Mean conductance as a function of the total flux 
(j) (a) and of the energy e. In (a), data are shown for e = 
and £ = 0.005. In (b), (g) is shown for (p = and = 8 for 
the RRH model (circles) and for the RF model at p = 0.3 (tri- 
angles). All simulations were done for A'' = 45 and L = 800. 
The average was performed over 10* realizations of th e dis- 
order. The dashed line in (b) is taken from Eq. (5.23) with 
/3 = 1 and cr^ = 1.5 x lO^e. 



Numerical results for the average conductance are 
shown in Fig. ||. All data shown are for the same length 
L ~ 800 and for the same number of channels = 45. 
For weak disorder one can ignore the and e dependence 
of the mean free path £, and hence of the Drude term in 
the conductance. The only effect of a variation of e or 
(j) is thus to change the symmetry of the quantum wire, 
which affects the weak-localization correction to the con- 



ductance. According to Eq. ( 5.25 ), we expect a nonzero 
weak-localization correction Sg in the standard orthog- 
onal symmetry class, i.e., for cj) — and £ ^ 0, while 
Sg = if time-reversal symmetry is broken ^ 1) or 
if the chiral symmetry is present (e = 0). This behav- 



ior is confirmed in Fig. |9[ However, qua ntita tively, the 
numerical results differ ^ 30% from Eq. ( [5^ ). In addi- 
tion. Fig. H shows a small e dependence of (g) at large 
magnetic fields that cannot be accounted for within the 
theory of Sec. In particular, note the cusp-likc struc- 
ture at small e in the = 8 data in Fig. This effect 
seems to be too large to be explained by a spurious e- 
dependence of the mean free path £. Since the feature 
at small e is suppressed at larger lengths L, a possible 
cause might be a contact resistance effect. (Contact re- 
sistance is known to play a role for disordered normal- 
metal-superconductor junctions, when the particle-hole 
degeneracv4s destroyed by a finite voltage or by a mag- 
netic field.C2l) As we discussed in the previous subsection, 
other causes for the discrepancy between theory and nu- 
merical simulations may be the fact that the disorder is 
not small, or that the parameter rj is not known. 



VII. CONCLUSIONS 

In the vicinity of the band center, the physics of lo- 
calization in a quantum wire with chiral symmetry ex- 
hibits differences with respect to the case of a quan- 
tum wire in one of the standard symmetry classes. The 
most prominent differences are observed in the localized 
regime. For wires with chiral symmetry (as is the case 
with off-diagonal disorder), at the band center, the statis- 
tics of the conductance depends sensitively on the parity 
in the number of transmission channels A'^. For odd N, 
the band center represents a critical point that is charac- 
terized by the absence of exponential localization. The 
logarithm of the conductance is not self-averaging and 
the mean conductance or its variance decay algebraically 
with the length L of the wire. For even N , exponential 
localization takes place with a self-averaging localization 
length ^ that does not depend on the presence or ab- 
sence of time-reversal and spin rotation symmetry. As 
the energy is tuned away from the band center, the sys- 
tem crosses over to the standard universality classes: The 
parity effect disappears and the localization length ac- 
quires a dependence on the presence or absence of time- 
reversal and spin rotation symmetry. In the presence 
of time-reversal symmetry, the localization length 1^ for 
even N is the same with and without chiral symmetry. 
Our numerical simulations indicate that the crossover is 
non-monotonous: in the crossover between the chiral and 
standard symmetry classes, ^ differs from the values in 
the pure symmetry classes. A complete theoretical de- 
scription of this crossover is still lacking. 

In the diffusive regime, the differences between the chi- 
ral and standard universality classes are less pronounced. 
They show up in quantum interference corrections to the 
classical (Drude) conductance, which is the same in both 
cases. We have found that in the chiral classes, weak- 
localization corrections to the mean conductance (g) and, 
more generally, to the density of transmission eigenval- 
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ues vanish at the band center. This is the quasi-one- 
dimensional counterpart of a similar observation made 
by Gade and Wegner in their study of two-dimensional 
disordered systems with chiral symmetry.Ell The conduc- 
tance fluctuations are twice as large at the band center 
relative to the standard universality classes, compatible 
with the presence of an extra symmetry in the system. 
We have calculated these quantum interference correc- 
tions as a function of energy, for the entire crossover from 
the chiral universality class at the band center e = to 
the standard unitary classes for e far away from 0. The 
theoretical predictions for this crossover agree qualita- 
tively with numerical simulations though there remains 
sizable deviations between theory and numerics of the 
order of 10-30%. 

While the chiral symmetry classes have received an 
enormous amount of theoretical attention (see the intro- 
duction of this paper for a brief summary), there are 
several hurdles to take before a chiral quantum wire can 
be realized in practice. Besides the effect of electron- 
electron interactions, which is not taken into account 
here, the main obstacle is the fact that the chiral sym- 
metry is very fragile, since it is easily broken by, e.g., 
on-site random energies, next-nearest-neighbor hopping, 
or a small shift of the chemical potential, which will drive 
the system away from the chiral symmetric band center. 
Our calculation of the quantum interference corrections 
in the crossover from the chiral symmetry classes to the 



standard ones can be seen as a first and necessary step 
to tackle the latter obstacle. 
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APPENDIX A: SCALING EQUATIONS 

In this appendix, we present the scaling equations 
needed to calculate the crossover for the weak localiza- 
tion corrections of the conductance, shot noise, and the 
universal conductance fluctuations of the conductance as 



done in Sec. V_C. The notation was defined in Eq. ( 5.11 ) 
and we use the short-hands 

R = i?io, i?2 = -Rioio, R3 = -Rioioio- 
Equations needed up to corrections of order N'^ are 



J 



{R) ^N{N^2R) + {Rf + ^ {N - R^o ~2R- R*o + R.ooo + ^tooo + ^2) , 
lldL{R2) = ^{R-R2)iN-{R)) 



(Al) 
(A2) 



2-/3 
/3 



(^00 + 4i? + Rqq 8i?2 — 4i?jgQo — 4i?Jooo + 4-^3 + 2i?ioiooo + ^^R^q^qoo + -^100100 ~t -"-looioo/ 1 



7^. 
J' 



— ^Lvar R- 4 ((i?) - N) var R + (R^„ + i?^^ + 2R- 2i?iooo ^ 2i?^ooo ^ 4i?2 + -Riooioo ^ ^"-looioo ^ ^^-^a/ 



2i?, 



2-/3 
/3 



(^00 + -^00 + 2i? — 2i?jooo ~ 2i?^Qoo — 4i?2 + -Riooioo + ^100100 + 2-R3) 



(A3) 



Equations needed up to corrections of order N are 



7^ 

-^Ol (i?oo) 



J 

J 



dL{R 



1000/ 



f3 



(i?oo) + A^ (A^ - 2i?oo> + (^oo> , (A4) 
(i?iooo) + N (2i?oo + 2R- 4i?iooo) - {R,, + 2R- 2R,,,,) {RJ - (i? - 2R,,J (R) , (A5) 



9l (^1100) - ^ {Rqo + 2i? + i?flo - 4-Riioo) - (2^^ - -Rnoo) (-^00) - (2^oo - 2^1100) (R) + (^1100) (^00) ' (A6) 



7£ 
J 



dL {R3) = 3 (2N - R3) + {R~ 4i?2 + 2i?3) (R) + (i?2 



(A7) 



(i?ioiooo>+2(A^(2i?i 



-R9 — 3R 



101000/ 



{R - i?iooo - R2 + Ri 



(Roo)) 
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